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Topological properties lie at the heart of many 
fascinating phenomena in solid state systems such 
as quantum Hall systems or Chern insulators. 
The topology can be captured by the distribution 
of Berry curvature, which describes the geometry 
of the eigenstates across the Brillouin zone. Em¬ 
ploying fermionic ultracold atoms in a hexagonal 
optical lattice, we generate topological bands us¬ 
ing resonant driving and show a full momentum- 
resolved measurement of the ensuing Berry cur¬ 
vature. Our results pave the way to explore 
intriguing phases of matter with interactions in 
topological band structures. 

Topology is a fundamental concept for our understand¬ 
ing of many fascinating systems that have recently at¬ 
tracted a lot of interest, such as topological supercon¬ 
ductors or topological insulators, which conduct only at 
their edges [T]. The topology of the bulk band is quanti¬ 
fied by the Berry curvature [2] and the integral over the 
full Brillouin zone is a topological invariant, called the 
Chern number. According to the bulk boundary corre¬ 
spondence principle, the Chern number determines the 
number of chiral conducting edge states [T]. While in a 
variety of lattice systems ranging from solid state sys¬ 
tems to photonic waveguides and even coupled mechani¬ 
cal pendula, edge states have been directly observed El- 
|7] , the underlying Berry curvature as the central measure 
of topology is not easily accessible. In recent years, ultra¬ 
cold atoms in optical lattices have emerged as a platform 
to study topological band structures and these sys¬ 
tems have seen considerable experimental and theoretical 
progress. Whereas in condensed matter systems, topo¬ 
logical properties arise due to external magnetic fields 
or intrinsic spin-orbit coupling of the material, they can 
in cold atom systems be engineered by periodic driving 
analogous to illuminated graphene [10] . Interestingly, the 
resulting Floquet system can have totally new topologi¬ 
cal properties m- The driving can, for example, be re¬ 
alized by lattice shaking HMD or Raman coupling nu¬ 
lls] with high precision control in a large parameter space. 
In particular, the driving can break time-reversal symme¬ 
try [T311IU1I6] and thus allows for engineering non-trivial 
topology jniiiH]. In quantum gas experiments, topolog- 



Figure 1. Engineering of the topology by band dress¬ 
ing. (a) Three laser beams intersecting under 120° interfere 
to form a tunable honeycomb lattice |26| with a variable off¬ 
set energy HAab between A and B sites, (b) Circular shak¬ 
ing of the lattice (shaded inset on the right) with frequency 
V resonantly {u ~ vab) couples the two lowest bands, which 
are plotted (red and blue lines) along a high symmetry path 
(KjMjK’jF: K)- (c) Dressed Floquet bands for different shak¬ 
ing amplitudes between 0 and 223 nm at a shaking frequency 
of i/ = 11 kHz. The dressed bands are calculated using Flo¬ 
quet theory |26| . (d) Sketch of the position of the topological 
defects in the undressed and dressed cases, illustrating the 
dramatic change of the topology. The arrow indicates the 
direction of the phase winding around the Dirac point, (e) 
2D dispersion relations of the bare and dressed bands show¬ 
ing a six-, respectively three-fold symmetry. In (c and e), 
the color represents the dressing (less color corresponds to 
stronger dressing). 


ical properties have been probed via the Hall drift of 
accelerated wave packets nans], via an interferometer 
in momentum space |20L 121] and in charge pump experi¬ 
ments |22| but so far, the underlying Berry curvature was 
not revealed. 
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Here, we demonstrate a method |231 IM| that, for the 
first time, allows for a direct measurement of the Berry 
curvature with full momentum resolution. We perform a 
full tomography of the Bloch states across the entire Bril- 
louin zone by observing the dynamics at each momentum 
point after a projection onto flat bands. The topological 
bands are engineered by resonant dressing of the two low¬ 
est bands of an artificial boron nitride lattice and feature 
a rich distribution of Berry curvature. From the Berry 
curvature, being the fundamental building block of topol¬ 
ogy, other relevant quantities such as the Berry phase or 
the Chern number can easily be obtained. 

Our system consists of ultracold fermionic atoms in a 
hexagonal optical lattice formed by three interfering 
laser beams. With an appropriate polarization |26| . a 
variable energy offset HAab between the A and B-sites 
(Fig. la), which breaks inversion symmetry, can be en¬ 
gineered. With the emerging band gap Hvab, the Dirac 
points at K and K’ become massive and for a large offset, 
the bands are flat (Fig. lb) [^. This is a key ingredient 
for our tomography, since the flat band acts as the refer¬ 
ence frame in which we reconstruct the eigenstates. Now, 
as a central experimental feature, we can dress the bands 
via circular shaking. In Fig. Ic, we show the dressed Flo- 
quet bands for near resonant circular shaking for different 
accessible driving amplitudes. Apart from the dramatic 
change in the dispersion relation, the topological prop¬ 
erties of the bands are changed. This manifests itself in 
the creation of a new Dirac point at the F-point and the 
annihilation of a Dirac point at the K point (Fig. Id). A 
three-fold symmetry also becomes visible in the disper¬ 
sion relation (Fig. le). 

The topological properties are not captured by the 
mere dispersion relation but by the Berry curvature, 
which describes the winding of the eigenstates across the 
Brillouin zone. Therefore, a complete tomography of the 
eigenstates of a Bloch band is mandatory for a measure¬ 
ment of the Berry curvature. The key idea behind our 
tomography is to reconstruct the eigenvectors from dy¬ 
namics after a projection onto flat bands [23 • Consider 
the Bloch sphere (Fig. 2a) whose poles are given by |k. A) 
and |k, 5), which are the Bloch states restricted to the A- 
respectively B-sublattice. The lower band can be writ¬ 
ten as |k) = sin(0k/2)|k. A) — cos(0k/2) exp(i0k)|k, H) 
and after a projection onto flat bands, the state oscil¬ 
lates around |k, H) with the frequency given by the 
energy difference of the flat bands. Then, the momentum 
distribution after time-of-flight is given by 


n(k, t) = /(k) [1 - sin(6lk) cos((()k + 27rr'kt)] (1) 


from which both 0k ^-nd 0k can be easily obtained, yield¬ 
ing the desired tomography of the eigenstates for each 
quasimomentum. Our method hence allows for a direct 


reconstruction of the Berry curvature according to 

r!(k) = ^ X • h (2) 

with h = (sin(0k) cos(0k), sin(0k) sin(0k), cos(0k)) |23]- 

Our experimental sequence for this state tomography is 


a 



b 
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Figure 2. Dynamical measurement of the topology of the 
dressed bands, (a) Illustration of the experimental protocol. 
We start in a band insulator in the lowest undressed band 

(1) and adiabatically ramp up the coupling strength (shaking 
amplitude) (ii). Now the Bloch state at each quasimomentum 
is a superposition of |k. A) and |k, B), pointing in different di¬ 
rections on the Bloch sphere. When the dressing is switched 
off (iii), the state is projected onto the flat bands and ro¬ 
tates on the Bloch sphere with the band difference u^. After 
time-of-flight, this yields an oscillation of the density at each 
momentum, (b) Experimental momentum distributions for 
different hold times after the projection onto flat bands. The 
color is the atomic density in momentum space (after time- 
of-flight), the hexagon marks the first Brillouin zone. Note 
that the clearly visible oscillations of the density at differ¬ 
ent momenta stem from interference of the two bands (Eq. 

(2) ). (c) Oscillations at different quasimomenta. The solid 
lines are sinusoidal fits. From the amplitude and phase of 
the oscillation, we reconstruct the dressed state according to 
Eq. (1). The experimental parameters are Vj = 15.15(15) Er, 
Aab = 11.45(11) kHz, 1 / = 11 kHz and a shaking amplitude 
of 223 nm j^S]. 
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Figure 3. Momentum-resolved measurement of the Berry cur¬ 
vature. (a) Amplitude (left column) and phase (central col¬ 
umn) obtained from the fits to the oscillations from Fig. 2, 
performed for each pixel in and around the first Brillouin zone 
(hexagon). From those fit results, we obtain, as our central re¬ 
sult, the momentum-resolved Berry curvature (right column), 
(see text for details). The experiment was performed for dif¬ 
ferent chiralities of the lattice shaking, destroying the Dirac 
point at either K or K’. (b) Theoretical results from a tight- 
binding Floquet calculation |26| using precisely the experi¬ 
mental parameters, yielding a very good agreement, (c) In 
the upper left corner, a zoom into a phase vortex (red square 
in (a)) is shown, illustrating the high momentum resolution 
of our method due to a pixel-wise evaluation of more than 
2800 pixels in the first Brillouin zone. The plots show the 
experimental Berry curvature (blue points) along the high- 
symmetry path in comparison with the theoretical (red solid 
lines) calculation. In (b) and (c), the Berry curvature is 
given in units of the inverse reciprocal lattice vector length 
|b| squared |26| . 


sketched in Fig. 2. We start with a cloud of 5 • 10^ 
single-component fermionic atoms forming a non¬ 
interacting band insulator in the undressed lattice. Due 


to a large offset between the A and B sites, IS.ab = 
11.45(11) kHz , the undressed bands are flat, such that for 
all quasimomenta |k) « |k, H). We adiabatically ramp 
up the shaking amplitude to 223 nm within 5 ms at a 
shaking frequency of 9 kHz and then ramp the frequency 
to 11 kHz within 2 ms. By suddenly switching off the 
dressing, we project onto the bare flat bands, such that 
Eq. ([^ can be applied. In Fig. 2b, we show typical 
time-of-flight images for different hold times in the flat 
bands. The images feature dynamics with very large 
contrast. The insets show time evolutions for different 
quasimomenta, revealing the pure sinusoidal oscillations 
with clearly distinct amplitudes sin(0k) and phases (/<k 
which are obtained by a simple fit to Eq. (Q . We observe 
very large and long-lived oscillations after the projection, 
yielding relative amplitudes of up to 0.8. Additionally, 
with more than 2800 pixels in the first Brillouin zone, 
the resolution in momentum space is very high. 

As the central result, we reconstruct the Berry curva¬ 
ture of the dressed band structure from these fits, which 
fully visualize the Bloch states as shown in Fig. 3a. The 
amplitude map features a pronounced three-fold symme¬ 
try, illustrating the breaking of equivalence between the 
K and K’ points. The amplitude has a maximum at the 
K point and is zero at the K’ and F points. Even more 
striking is the very distinct threefold symmetry of the 
phase map with nearly discrete values of 7r/3, tt, 57r/3. 
Where the amplitudes are zero at the K’ and F points, the 
phase map correspondingly displays vortices. The phase 
vortices are clear signatures of Dirac points which con¬ 
stitute topological defects. Furthermore, the data clearly 
shows that we annihilate the Dirac point of the undressed 
hexagonal lattice at K and created a Dirac point in the 
dressed system at the F point, thereby strikingly chang¬ 
ing the topology of the band. The resulting Berry curva¬ 
ture is localized at the new Dirac points and also shows 
this clear three-fold symmetry. It has opposite sign at 
the two Dirac points which results from the opposite chi¬ 
rality of the phase vortices. By inverting the chirality of 
the shaking, we instead annihilate the Dirac point at the 
K’ point and invert the chirality of the phase windings, 
which also results in an inverted symmetry in the Berry 
curvature. All quantities agree well with a Floquet theory 
calculation (see Fig. 3b), based on a tight-binding model 
as described in m- In Fig. 3c, we plot the Berry cur¬ 
vature pixel-wise evaluated along a high-symmetry path, 
illustrating the very good agreement with the theory. 

As mentioned before, with the fully momentum- 
resolved Berry curvature, we can easily obtain the further 
relevant quantities like the Berry phase or Chern number. 
For a discussion of the respective Berry phases, see |26| . 
The integral over the closed area of the full first Brillouin 
zone must be quantized to 27r times the integer Chern 
number C. From our data we obtain C = 0.010(13) 
and C = —0.019(13) for the two different shaking chiral¬ 
ities |26| . which clearly confirms this quantization within 
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the experimental errors. Our data constitutes the first 
momentum-resolved measurement of the Berry curvature 
of a topological band structure, thereby providing a deep 
insight into the fascinating physics of topological sys¬ 
tems. Our measurements demonstrate that even when 
the global topology has Chern number zero, the distri¬ 
bution of Berry curvature can be very rich. 

Our measurement scheme can be readily extended to 
characterize bands with Chern numbers different from 
zero IHIIH]. In principle, one could start in a shallow lat¬ 
tice, where reaching non-zero Chern numbers is possible, 
and for the tomography project onto flat bands, which 
can be reached, e.g., by dynamical control over the offset. 
Our method for generating the topological bands is spin- 
independent and does not couple different spin states. It 
therefore can be extended to high-spin systems m or to 
strongly-interacting spin mixtures, which are expected to 
lead to interesting many-body phases |2 81150] . 

We acknowledge stimulating discussions with Andre 
Eckardt, Christoph Strater and Maciej Lewenstein. This 
work has been supported by the excellence cluster "The 
Hamburg Centre for Ultrafast Imaging - Structure, Dy¬ 
namics and Control of Matter at the Atomic Scale" and 
the GrK 1355 of the Deutsche Forschungsgemeinschaft. 
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A Tunable driven hexagonal lattice 

For the experiments presented in the main text, we 
employ a new lattice setup allowing for a variety of lat¬ 
tice geometries by controlling the polarizations of the in¬ 
terfering laser beams, following a proposal as described 
in m- The lattice potential is given by the sum of a 
honeycomb lattice Vp(r) generated by the p-polarization 
components of the laser beams and a triangular lattice 
14(r) due to the s-polarization components. By tuning 
the populations 9 and relative phases a between the s- 
and p-components, the relative depth and position of the 
triangular lattice with respect to the honeycomb lattice 
can be set. In particular, we can continuously transform 
our pure and totally symmetric honeycomb lattice into 
a boron-nitride lattice in which the A- and B-sites are 
no longer energetically degenerate and inversion symme¬ 
try is broken (see Fig. la in the main text and Fig. |^, 
as done for the experiments presented in the main text. 
The periodic driving of the lattice is realized by modulat¬ 
ing the frequencies of beams 2 and 3 with acousto-optic 
modulators which is equivalent to modulating the rela¬ 
tive phases with respect to beam 1, ip 2 and 1 ^ 3 . This 
allows for moving the lattice in real space on nearly ar¬ 
bitrary trajectories. In Fig. our tunable lattice setup 
is sketched. 


sin(03) exp(ia3) 6^+ 005(03) sinlOj) 6xp(ia3) 6^+ cos(0j) 



Figure 4. A tunable hexagonal lattice setup. The lattice 
is formed by three interfering laser beams intersecting under 
an angle of 120°. The periodic driving is realized by modu¬ 
lating the relative phases of the beam. The p- (in-plane) 
polarization interference pattern is a honeycomb lattice while 
the s- (out-of-plane) polarization interference pattern is in¬ 
verted and forms a triangular lattice. The relative phase a 
between s- and p-polarization shifts the triangular lattice with 
respect to the honeycomb lattice. The population 9 deter¬ 
mines the relative depth of the two lattices. For 9 = 7r/20 for 
all beams and 01 , 2,3 = 1/3(0, 2?!, dir), a boron-nitride lattice 
with distinct A and B sites is realized. The resulting disper¬ 
sion relation is plotted for the experimental parameters as in 
Fig. 2 of the main text, showing that the lower two bands 
(red, blue) are well separated from the other bands (gray). 


The lattice depths Vi stated in the main text are 
specified in units of the recoil energy Ej. = h^k\l2m 
with kL = 27r/A and the lattice laser wavelength is 
A = 1064 nm, Er/h « 4410 Hz. The complete lattice 
potential can be written as F(r) = V),(r) -|- I 4 (r) with 
1 ^ (r) = ViCOS^( 0 )(cos(bi-r-v 32 )+cos(b 2 -r-|-(( 92 -<i 33 ) + 
cos((bi -bb 2 ) -r —(/J 3 )) and 14(r) = —2Vi sin^(0)(cos(bi) • 
r -I- ai - 02 - 1^2 ) + cos(b 2 • r -I- 02 - 03 -I- (^2 - )+ cos ((bi -I- 
b 2 )-r-|-oi— 03 —(/J 3 )). Here, bi = ki—k 2 and b 2 = k 2 —k 3 
are the reciprocal lattice vectors defined via the wave vec¬ 
tors ki = fci,(0,1,0), k 2/3 = fci(±v^/2,-1/2,0). For 
the experiments presented in the main text we choose 
Vi = 15.15(15) Er, 9 = 7 r /20 and a, = {j - 1) • 27r/3. 
In Fig. SI, the resulting dispersion relation of the six 
lowest bands for the experimental parameters is plotted 
(lower right corner). Because the bandgap (11680 Hz) 
is large compared to the bandwidths of the two lowest 
bands (335 Hz) and (48 Hz), the bands can be consid¬ 
ered flat. The lowest two bands are energetically well 
separated from the other bands. The lattice depth is 
obtained from the fitted oscillation frequencies (Eq. 
1 of the main text) and the dominating systematic er¬ 
ror of 1% stems from the lattice depth inhomogeneity 
across the sample. In the transverse direction perpen¬ 
dicular to the plane spanned by the lattice beams, the 
potential is harmonic, thus forming an overall lattice of 
tubes. In addition to the lattice potential, a crossed 
dipole trap operated at a wavelength of 825 nm is su¬ 
perimposed, formed by a round beam perpendicular to 
the lattice direction with a 1 /e^-intensity beam radius 
of Wx^y = (70, 240) ^m and an elliptical beam propagat¬ 
ing in the x-direction with Wy^z = (70, 240) ^m. The 
measured trapping frequencies for the experimental pa¬ 
rameters are Vx,y,z = (83(4), 108(5),93(4)) Hz. 


Tight-binding description of the lattice Hamiltonian 


Since the two lowest bands are well separated from 
the other bands (see Fig. |^, a two-band tight-binding 
model jS2] can be employed to describe our system. The 
tight-binding Hamiltonian can be written as 
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H(k) 


f + E* cos(k • ai 
V Ei2tABexp(+zk-di) 


with the offset Aab between A- and B-sites and the 
next-nearest-neighbor hopping amplitudes tAA, ^bb on 
the diagonal and on the off-diagonal the nearest-neighbor 
hopping amplitude Iab coupling the A- and B-sites (see 

Fig.|^. 

The vectors connect neighboring A(B)-sites with 
each other while the vectors connect A- with neighbor¬ 
ing B-sites. For our perfectly symmetric hexagonal lat¬ 
tice, ai = — 47r/3A:L(0, 1 , 0), ai = —47r/3fcL(-\/3/2,1/2, 0) 
and as = a 2 — ai are the Bravais vectors, di = l/3(ai - 1 - 
^ 2 ), d 2 = l/3(—2ai -I- a 2 ) and ds = l/3(ai — 2a2). We 


) Ei 2 iAB exp(-ik • di) \ 
+ Ei 2 tBB cos(k • a;)/ 


obtain the tight-binding parameters Aab = 11-45 kHz, 
tAB = 657 Hz, tAA = 32 Hz and tsB = 0 Hz from a fit 
to the exact band structure shown in Fig. The asym¬ 
metry between the lower and upper band is captured by 
the differing next-nearest-neighbor hopping amplitudes 
tAA 7 ^ tSB- The uncertainty of Aab stated in the main 
text arises due to the systematic error of 1 % in the lattice 
depth (see the previous section). The lower band eigen¬ 
vectors |k) = sin(0k/2)|k. A) — cos(0k/2) exp(i(/k)|k, H) 
can be expressed in terms of the angles 0 k and (/k which 
are calculated according to: 


sin( 0 ) = 


|X)iUBexp(fk-d.,)| AAB/2-f Ej(^AA-<BB)cos(k-aj) 

-jCOS^Wj = - 


R 


R 


R = 


- arg tAB exp(ik • d^) j 

^Aab/ 2 -f '^{tAA - tBB) cos(k • ai)^ -t- tAB sin(k • d^) 


( 4 ) 



1/2 


Floquet description of the dressed system 


A periodically driven system can in Floquet theory 
be described by an effective (time-independent) Hamilto¬ 
nian Heff whose properties characterize the dressed sys¬ 
tem. The effective Hamiltonian is obtained from the 
unitary time-evolution operator U (T, 0) over one pe¬ 
riod T of the driving |9l |33], ^ log(U(T, 0)) = Heg. 
The time evolution operator is given by U (T, 0) = 
exp ^ jJ" H(t)dt^ and can be numerically obtained 
by replacing the integral with the time-ordered product 
U(r, 0) « U(T - nAt, T-{n+ l)At). Using 

the approximation U(t-|-At, t) « exp ^H(t) • At^ , the 
product can be calculated. We ensured convergence nu¬ 
merically by checking that increasing the amount of time 
steps does not change the results. Using the tight-binding 
description of the lattice Hamiltonian as described above 
(Eq. allows for a computationally cost-effective cal¬ 
culation of the effective, time-independent Hamiltonian 
whose eigenvalues describe the dispersion relation of the 
dressed bands. Using Eq. ( 0 . the angles describing the 
dressed (effective) eigenvectors on the Bloch sphere are 


obtained, which allows for a direct comparison with the 
experimental data as presented in Fig. 3 in the main 
text. 


Derivation of the data for the calculation of the 
Berry curvature 

To reconstruct the Berry curvature, we fit n(k, t) = 
/(k) [1 — Aexp(—t/r) cos( 0 k + 27rr'kt)] to each pixel, ac¬ 
cording to Eq. (1) in the main text. The fit parameters 
are the Wannier envelope /(k), the relative amplitude A, 
the damping constant r, the phase (/k and the frequency 
zzk which is given by the energy difference between the 
lower and upper undressed bands. The damping time 
typically is larger than five oscillation periods. We scale 
the experimentally measured amplitude by a factor of 
1.15, such that the maximum amplitude is 0.92, as ex¬ 
pected from theory. We attribute this scaling factor to 
two effects. First, both the harmonic trap coupling differ¬ 
ent quasimomenta and the inhomogeneous lattice depth 
lead to a dephasing and hence reduction of the amplitude 
of the oscillations. Second, due to unavoidable multi¬ 
photon excitations into higher bands which lead to par- 
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Figure 5. Sketch of lattice geometry. Arrows indicate 
hopping possibilities Iaa, tsB and Iab- 


tide loss during the ramp, we choose the ramping time 
to be not perfectly adiabatic. We quantified the latter 
by ramping back into the undressed system and a sub¬ 
sequent band mapping and found a relative population 
in higher bands of 13 %. This explains why the exper¬ 
imental amplitude is approximately 15 % smaller than 
expected from theory and justifies the rescaling. While 
the phase directly corresponds to the desired tight- 
binding parameter, the other parameter is obtained via 
61k = arcsin(d). When both the north and south hemi¬ 
spheres of the Bloch sphere are covered (as for non-zero 
Chern numbers), one can obtain the sign from an addi¬ 
tional band mapping measurement or by a careful anal¬ 
ysis of the zeros in the amplitude . Before fitting, we 
smooth n(k, t) along the time direction with a moving 
average of span 3. For the calculation of the Berry cur¬ 
vature according to Eq. (2), both 0k and and their 
derivatives are needed. While 0k and ^k are evaluated 
pixel-wise, their derivatives are obtained from Savitzky- 
Golay filtered (3''^ order polynomial and window size 15 
pixels) finite differences. The pixel-size is small compared 
to the length (|b| = 57(2) px) of the reciprocal lattice 
vectors, allowing for a very high resolution in momentum 
space. As shown in Fig. 3 in the main text, the exper¬ 
imentally reconstructed Berry curvature is in very good 
quantitative agreement with the theoretical calculation. 
As stated in the main text, the Berry phase and Chern 
number can be directly obtained from the Berry curva¬ 
ture. The Chern number which is given by the integral 
of the Berry curvature over the first Brillouin zone is in 
our case simply the sum of Berry curvature over all pixels 
within the first Brillouin zone (multiplied by the area of 
the Brillouin zone and divided by 27r). The errors stated 
for the Chern numbers result from the uncertainty in the 
length of the lattice vector and an uncertainty of 2 px in 


the position of the first Brillouin zone. Since the Berry 
phase around a loop is the integral of the Berry curvature 
over a closed area, it simply is the sum of Berry curva¬ 
ture over the pixels in that area multiplied by the area. 
Hence, the Berry phase can be directly obtained from our 
data for arbitrary loops, proving the powerfulness of our 
method. 
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